Atomistic Modeling of Quaternary Alloys: Ti and Cu in NiAl 
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The change in site preference in NiAl(Ti,Cu) alloys with concentration is examined experimentally 
via ALCHEMI and theoretically using the Bozzolo-Ferrante-Smith (BFS) method for alloys. Results 
for the site occupancy of Ti and Cu additions as a function of concentration are determined experimen- 
tally for five alloys. These results are reproduced with large-scale BFS-based Monte Carlo atomistic 
simulations. The original set of five alloys is extended to 25 concentrations, which are modeled by 
means of the BFS method for alloys, showing in more detail the compositional range over which 
major changes in behavior occur. A simple but powerful approach based on the definition of atomic 
local environments also is introduced to describe energetically the interactions between the various 
elements and therefore to explain the observed behavior. 


I. INTRODUCTION 

THE influence of ternary alloying additions on the struc- 
ture of binary ordered alloys is an area in which theory and 
experiment have kept pace with each other. Experimental 
techniques provide accurate and detailed data that are rather 
easily matched by current theoretical methods. Abundant 
work in this area has resulted in a thorough understanding 
of site occupancy and the underlying phenomena affecting 
the structure of materials. In spite of this promising starting 
point, severe limitations are almost immediately found if 
similar information is needed from more complex systems. 
That is definitely the case for the analysis of several minority 
elements in an otherwise simple (i.e„ binary) ordered alloy. 
It is only recently that some progress has been made in such 
analysis, at a fundamental level, due to the emergence of 
more accurate experimental techniques or, in the case of 
theoretical work, more computer power to attack complex 
problems.* 1 ” 43 Quantum approximate methods offer a new 
tool for investigating such problems by introducing simple 
formulations that translate into computationally efficient 
calculations. 

In this article, we move one step forward, both theoreti- 
cally and experimentally, in understanding the role of simul- 
taneous alloying additions by studying a quaternary system, 
a NiAl-based alloy with Ti and Cu alloying additions. It is 
obvious that as the number of dements increases, the amount 
of available experimental evidence decreases, thus assigning 
modeling techniques a growing role in filling the gaps 
in our knowledge. The purpose of this work is to highlight 
this point, by presenting recent experimental work, using 
Atom Location by Channelling Enhanced Microanalysis 
(ALCHEMI/ 53 on NiAl(Ti,Cu) alloys.* 63 The experimental 
data were used to verify the modeling results and to validate 
the quantum approximate method used, the Bozzolo- 
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Ferrante-Smith (BFS) method for alloys. [7J Once the model- 
ing approach was validated, it was extended to a much 
broader set of compositions. We conclude with an analysis 
of the energetics of the system for a wide range of concentra- 
tions, including those studied experimentally. In doing so, 
we attempt to show how powerful modeling techniques sup- 
plement necessarily limited experimental data while at the 
same time providing global explanations of the observed 
behavior. To this end, this article is organized as follows. 
Section II describes the A.LCHEMI results for five B2-struc- 
tured NiAl(Ti,Cu) alloys (Ni 5 oAl ( 47 _ x /Ti 3 Cu*, with X = 1, 
3, and 6, and Ni^o-^ALyT^Cu*, with X = 1 and 3), as 
shown in Figure 1. The site preference behavior of the Ti 
and Cu alloying additions as a function of Ni and A1 concen- 
tration was determined. Section III describes the modeling 
effort, starting with the presentation of the results of Monte 
Carlo-Metropolis large scale atomistic simulations of those 
five experimental alloys, thus validating the use of the BFS 
method for alloys in this particular system. The modeling 
effort is then expanded to include detailed analytical BFS 
calculations for a wider range of concentrations, in order to 
provide much needed detail on the energetics and the 
resulting structure of the ground states. This is followed, in 
Section IV, with the introduction of a simple approach based 
on the BFS analysis for describing the energetics of the local 
atomic environments of the substitutional defects. Conclu- 
sions are drawn in Section V. 

IL ALCHEMI ANALYSIS OF NiAl(TI ? Cu) 
ALLOYS 

A. Experimental Procedure 

The alloys in this study were created as variations of the 
base composition Ni 5 oAl 4 7 Ti 3 (at. pet), as shown in Figure 
i . In alloys 1 through 3, the amount of A1 replaced by Cu 
was 1, 3, and 6 at. pet, respectively. Alloys 4 and 5 contain 
1 and 3 at. pet Cu, respectively, added to replace Ni. The 
alloys were given a heat treatment consisting of a final 
homogenization treatment of 32 hours at 1644 K and aged 
for 6 hours at 1255 K, followed by air cooling. Transmission 
electron microscope (TEM) samples were prepared from 
3-mm disks mechanically ground to 250 jam thickness. The 
ground disks were electropolished using a solution of 70 pet 
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Fig. 1 — Composition of the NiAl-Ti-Cu alloys modeled. Those for which 
experimental data exist [5] are denoted with filled circles. The numbers inside 
the circles indicate the concentration of Cu (in at. pet). The horizontal and 
vertical axes indicate the concentration of A1 and Ni, respectively. The Ti 
concentration was held constant at 3 pet. The solid circles denote the alloys 
studied analytically with BFS. 


ethanol, 14 pet distilled water, 10 pet butylcellosolve, and 
6 pet perchloric acid at 20 V and 258 K. Energy-dispersive 
X-ray spectroscopy (EDS) was performed on a JEOL* 

*JEOL is a trademark of Japan Electron Optics Ltd., Tokyo. 

201 OF TEM operating at 200 kV. The X-ray spectra were 
acquired and analyzed using a NIST Desktop Spectrum Ana- 
lyzer (DTSA version 2.5.1). A cold stage specimen holder 
was used to minimize thermal vibrations, which may reduce 
the electron channeling required for the analysis. An area 
approximately 200 nm in diameter and 150-nm thick was 
used for the analyses with a beam convergence angle of 
3.5 mrad. 


B. Electron Channelling-Enhanced Microanalysis 

The site occupancies were determined using the TEM- 
EDS technique of ALCHEMI. [51 This technique requires 
the acquisition of EDS spectra under carefully controlled 
dynamical electron diffracting conditions in a crystal. In the 
case of B2-ordered NiAl, an optimum condition occurs in 
a two-beam diffracting condition near g — (100): [5} By alter- 
ing the diffraction conditions, the electron intensity along 
compositionally distinct planes can be controlled. This effect 
can be observed by acquiring EDS spectra at several diffract- 
ing conditions and measuring the relative changes in X-ray 
intensity for Ni and Al. If Ti or Cu is situated predominantly 
on one type of lattice site (either Ni or Al), its X-ray intensity 
will be altered in a like manner to the host element as 
diffracting conditions change. To explicitly calculate site 
occupancies, the diffraction-induced changes in generated 
X-ray intensities must be directly connected between 
alloying elements and host elements. Several methods of 
calculating site occupancies have been published in the liter- 
atureJ 8 ' 121 In this analysis, the equations in Reference 13, 
described subsequently, were used to calculate the percent- 
age of Ti and Cu on Al and Ni sites. In the case of Ti, X- 
ray generation is related to Ni and Al intensity by 

Wi = ^Ni^Ni + a AI^Al + CVi [1] 


where I x is the measured X-ray intensity of element x. The 
coefficients a Ah a Ni , and C Tl can be fit using multiple regres- 
sion [15] by collecting a set of l Tv / Ni , and / A1 values for several 
diffracting conditions. Once the coefficients in Eq. [1] have 
been determined, the relative amount of Ti on Ni and Al 
sites can be calculated. Following the procedure outlined 
in Reference 13, which assumes that dynamical correction 
factors are independent of diffracting conditions, and assum- 
ing that the only relevant dynamical correction is for Al-K, 
the relative amount of Ti on Ni sites, / Ni , is given by 
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By accounting for the effect of Al-K delocalization in this 
analysis, the results are slightly different than those presented 
previously, 1 [6] where the second term is used to compensate 
for delocalization of generated X-rays. These values can be 
normalized under the assumption that Ti can only reside on 
Ni or Al atomic sites. 

Errors in this method of site occupancy measurement can 
be attributed to several factors. Anti-site defects, where Ni 
atoms do not exclusively reside on Ni lattice sites, will 
directly affect the measured channeling effect. Delocaliza- 
tion, where characteristic X-ray generation is not completely 
limited to lattice points, reduces the channeling effect for 
lower energy X-rays. fl4] Standard sources of errors associ- 
ated with EDS microanalysis in the TEM must also be 
considered. The error in the site occupancy calculation for 
this analysis is limited to the variables in each equation. 
Most of the errors in the calculations were due to the uncer- 
tainties in the multiple regression coefficients; however, the 
statistical errors in the measured X-ray intensities were also 
included. Assuming all variables in Eqs. [2] and [3] are 
independent, the uncertainty in the site occupancies can be 
estimated by the following equations: 
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C. Experimental Results 

After collecting the EDS spectra from several diffracting 
conditions (including kinematic) for each alloy, the spectra 
were background subtracted and the intensities were meas- 
ured for the Al-K, Ni K*, Ti K„, and Cu K* peaks. To 
illustrate the channeling effect, two EDS spectra collected 
at different diffracting conditions from alloy 3 are displayed 
in Figure 2. The spectra are normalized with respect to the 
Ni K a peak. It is clear that the Al-K, Ti K a and Cu K a peaks 
are dramatically changed relative to the Ni K a peaks. This 
initial examination of spectra indicates that the channeling 
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Fig. 2 — EDS spectra from Ni 50 Al 4 iTi 3 Cu 6 alloy at two strongly diffracting 
conditions near g = (100). 


Table I. Percentage of Ti and Cu on A1 Sites 


Alloy 

Ti 

Cu 

NisoAl^T^Cui 

93.1 ± 5.2 

77.1 ± 6.4 

Ni 50 Al4 5 Ti 3 Cu2 

90.1 ± 10.0 

91.8 ± 17.9 

Ni 5 oAl44Ti 3 Cu3 

95.5 ± 8.0 

88.1 ± 13.1 

Ni49Al 47 Ti 3 Cu 1 

85.9 ± 6.4 

29.6 ± 8.3 

Ni47Al 47 Ti 3 Cu 3 

87.3 ± 15.9 

6.8 ± 18.7 


effect is strong, and both Ti and Cu preferentially reside on 
A1 lattice sites in this alloy. 

Similar analysis was performed on all five alloys and the 
site occupancies for Ti and Cu (Eqs. [2] and [3]) are listed 
in Table I, along with the calculated uncertainties. The results 
demonstrate that the Cu content on A1 sites is strongly 
dependent upon the stoichiometry of the alloy, whereas the 
Ti strongly prefers A1 sites for all alloys. A major portion 
of the uncertainties listed in Table I was due to statistical 
error in the C Ti constant, which was included to incorporate 
a delocalization correction for the Al-K X-ray line (Eq. 
[3]). It also introduced a significant amount of error in the 
calculation. Nevertheless, the calculation indicates that Cu 
substitutes for either Ni or A1 sites, depending on 
stoichiometry. 


III. THEORY 

A. The BFS Method for Alloys 

The BFS method for alloys has been extensively used for 
the analysis of ordered intermetallic alloys due to its ability 
to describe multicomponent systems with the same simplic- 
ity and accuracy as binary systems. [7, 16-181 Unlike other 
quantum approximate methods of similar computational 
characteristics, BFS benefits from its universal parameteriza- 
tion, which allows it to deal with NiAl, NiAlTi, NiAlCu, 
and NiAlTiCu alloys with the same set of parameters, regard- 
less of concentration or composition, as long as we are 
dealing with a bcc-type structure, including B2, Heusler, 
and other bcc derivative compounds. 


In this section, we provide a brief description of the opera- 
tional equations of BFS. The reader is encouraged to seek 
further details in previous articles where a detailed presenta- 
tion of the foundation of the method, its basis in perturbation 
theory, and a discussion of the approximations can be 
found. 171 The BFS method provides a simple algorithm for 
the calculation of the energy of formation of an arbitrary 
alloy (i.e., the difference between the energy of the alloy 
and that of its individual constituents). In BFS, the energy 
of formation is written as the superposition of elemental 
contributions of all the atoms in the alloy: 

AH = y 8, [6] 

i 


For each atom, we partition the energy s, into two parts: 
a strain energy and a chemical energy contribution. The BFS 
strain energy relates only to the atomic positions of the 
neighboring atoms to atom i, regardless of their chemical 
identity. For its calculation, we use the actual geometrical 
distribution of the atoms surrounding atom i, computed as 
if all of its neighbors were of the same species as atom i. 
Its calculation is then straightforward, even amenable to 
first-principles techniques. In our work, we use equivalent 
crystal theory (ECT) [191 for its computation, due to its proven 
ability to provide accurate and computationally economical 
answers to most general situations. 

The chemical environment of atom i is considered in the 
BFS chemical energy contribution, where the surrounding 
atoms maintain their identity but are forced to occupy equi- 
librium lattice sites corresponding to the reference atom i. 
Building on the concepts of ECT, a straightforward approach 
for the calculation of the chemical energy is defined, properly 
parameterizing the interaction between dissimilar atoms. [7] 

Thus defined, the BFS strain and chemical energy contri- 
butions take into account different effects, Le., geometry 
and composition, computing them as isolated behaviors. A 
coupling function, g,-, restores the relationship between the 
two terms, defined in such a way as to properly consider the 
asymptotic behavior of the chemical energy, where chemical 
effects are negligible for large separations between dissimilar 
atoms. The contribution by atom i to the energy of formation 
is then e* = sf + g^f. The strain energy contribution, sf, 
is obtained by solving the ECT perturbation equation: 


NRf l e ~ a ' R i + MR^eV^P 


= 2 rpe' (ai+S(r ^ [7] 


where N and M are the number of nearest-neighbors (NN) 
and next-nearest neighbors (NNN), respectively, and where 
p , /, a , and A are ECT parameters that describe element i 
(Reference 19 provides definitions and details), r denotes 
the distance between the reference atom and its neighbors, 
S(r) describes a screening function, [I91 and the sum runs over 
NN and NNN. This equation determines the lattice parameter 
of a perfect equivalent crystal where the reference atom i 
has the same energy as it has in the geometrical environment 
of the alloy under study. The terms R j and R 2 denote the 
NN and NNN distances in this equivalent crystal. Once the 
lattice parameter of the (strain) equivalent crystal, a s , is 
determined, the BFS strain energy contribution is computed 
using the universal binding energy relation of Rose et al} 20] 
ef = Ec(l “ (1 + af*)e ~ a * ), where E' c is the cohesive 
energy of atom i and where the scaled lattice parameter 
af* is given by af* = (<?// ( )(af* ~ af e ), where q is the ratio 
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between the equilibrium Wigner-Seitz radius and the equi- 
librium lattice parameter a l e . 

The BFS chemical energy, ef, is obtained by a similar 
procedure. The surrounding atoms retain their chemical iden- 
tity, but are forced to be in equilibrium lattice sites of an 
otherwise monatomic crystal i. The BFS equation for the 
chemical energy is given by 


NR^e~ a ^ + MR Per*)* 2 


[ 8 ] 


= I [N ik rfe' aikrk + M^eT^ 

k 


where N ik and M ik are the number of NN arid NNN of species 
k of atom i. The chemical environment surrounding atom i 
is reflected in the parameters a ik , given by a ik = a t + A*,-, 
where the BFS parameters A (a perturbation on the single- 
element ECT parameter a,-) describe the changes of the wave 
function in the overlap region between atoms i and k. Once 
Eq. [8] is solved for the equivalent chemical lattice 
parameter ap, the BFS chemical energy is then ep = 
y i E i c (\ - (1 4- af)e~ a ^ ), where y = 1 if af* > 0 and y = 
- 1 if af* < 0. The scaled chemical lattice parameter is 
given by ap* = (q/li)(af* - a*,). Finally, as mentioned pre- 
viously, the BFS chemical and strain energy contributions 
are linked by a coupling function g h which describes the 
influence of the geometrical distribution of the surrounding 
atoms in relation to the chemical effects and is given by 
gi ~ exp The pure element parameters a e , E c , /, a , 

and A, computed using the linear-muffin-tin orbital method 
(LMTO) [211 in the atomic sphere approximation, and the BFS 
parameters, A ab and A ba , for A and B = Ni, Al, Ti, and 
Cu, used in this study are listed in Reference 17. These same 
parameters are used in any calculation involving these four 
elements, regardless of the type of defects present (including 
surfaces) and the composition of the alloy, as long as the 
final crystal structure remains bcc. 


B. BFS Modeling 

The first step in our modeling effort consists of testing 
the ability of the theoretical method to reproduce the experi- 
mental results. In this case, it is necessary to perform com- 
puter experiments for the same five concentrations studied 
via ALCHEMI, as shown in Figure 1. Clearly, no single 
theoretical technique can reproduce all the subtleties of an 
experiment, much less simultaneously deal with ail the 
length and time scales involved. It is not strictly necessary, 
however, to perform a perfect theoretical simulation of the 
experiment. For the purpose of this work, it suffices to 
implement a ground state search where a large collection of 
atoms is allowed to evolve to its lowest energy state follow- 
ing a predetermined temperature-dependent process. This is 
best achieved by means of Monte Carlo -Metropolis simula- 
tions/ 161 Due to the number of elements involved in these 
calculations, a full Monte Carlo treatment should include 
individual atomic relaxations and the ability of the system 
to evolve to a multiphase structure. However, previous work 
on the determination of the site occupancy of additions to 
NiAl fl6 ’ 17,22J indicates that, for the range of concentrations 
studied, it suffices to consider only global relaxations, where 
only the lattice parameter of the alloy (/.£., computational 
cell) is varied until the energy is minimized, thus introducing 


a substantial simplification in the calculations. It should be 
noted that in the absence of a satisfactory foundation for 
this approximation, a full treatment of relaxations must be 
performed. 

L. Large scale simulations 

Monte Carlo-Metropolis simulations using the BFS 
method are performed in a sufficiently large atom population 
(1024 atoms). In these simulations, pairs of atoms at NN 
distance are allowed to switch places with a probability exp 
(— AEfkT), where k is the Boltzmann’s constant, T is the 
temperature, and A E is the difference in energy of the cell 
before and after the switch. The changes in atomic distribu- 
tion are allowed to continue until the total energy of the cell 
stabilizes. Every computational cell has been subject to the 
same temperature treatment, where an initial random distri- 
bution of atoms is steadily cooled from over 2000 K to 
room temperature. 

2. Verification of Monte Carlo results with experiment 

The results for the five alloys that were studied experimen- 
tally^ 1 are shown in Figure 3. Besides the obvious insight 
that can be obtained from direct visual inspection, valuable 
information can also be extracted from the coordination 
matrix of the final, stable state for each cell, where the 
matrix element a mn represents the probability that an atom 
m has an atom n as a NN. 

Table II displays the coordination matrices for the five 
alloy concentrations studied. As long as the cell follows the 
basic B2 ordering (in this case, Ni and Al atoms occupying 
their own sublattices), then <z mNi and a mAl can be taken as 
an approximate measure of the likelihood that an atom m 
occupies a site in the Al or Ni sublattice, respectively. If 
P[m(n)] denotes the probability of an atom m occupying a 
site in the n sublattice, then the previous statement can be 
written as a mH] ~ P[w(Al)] and a mM ~ P[m(Ni)]. Large 
values of the diagonal elements in each matrix can indicate 
either antiphase boundaries or the presence of precipitates 
(Le., a large value of n Ni Ni indicates that many Ni atoms are 
at NN distance, which would be highly unlikely in a NiAl 
alloy, where, if perfectly ordered, a NiNi = 0). If the diagonal 
elements are small, then the off-diagonal elements can be 
taken as a good approximation of the site preference. In 
other words, if a NiNi and n A1A1 are small, the closer the cell 
is to a highly ordered state, which translates into a TlNi being 
a true measure of the likelihood of finding Ti in an Al site. 

In agreement with experimental results and the conclu- 
sions reached in Reference 6, the simulations show that Cu 
does change site preference in NiAl, depending on the ratio 
of Ni to Al, as seen in the variations of a CuN j in Table II. 
For the five alloys shown, P[Ti(Al)] > P[Ti(Ni)], whereas 
P[Cu(Al)] > P[Cu(Ni)] for Ni-rich alloys and P[Cu(Ni)J > 
P[Cu(Al)] otherwise. Figure 3(f) shows a comparison of 
the computed site occupation probabilities PfTi(Al)] and 
P[Cu(Al)] with those determined experimentally/ 61 clearly 
showing the agreement between experiment and theory and 
providing a great degree of confidence in the BFS-based 
simulations. 

The possibility of a general trend in site occupancy as a 
function of NiAl concentration warrants additional simula- 
tions for a range of concentrations that includes the five 
studied experimentally (Figure 1). Figure 4 shows P[Ti(Al)] 
(Figure 4(a)) and P[Cu(Al)] (Figure 4(b)) as a function of 
Al, Ni, and Cu concentration. As expected, P[Ti(Al)] is 
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Fig. 3 —(a) through ( e ) Final structures of the Monte Carlo/Metropolis/BFS simulations of the five alloys studied in Ref. 6. White, light gray, dark gray, 
and black spheres denote Ni, Al, Cu, and Ti atoms, respectively. The 1024 atom cells are stretched along the vertical axis to guide the eye. (/) Comparison 
of the occupation probabilities P[Ti(Al)] and P[Cu(Al)] from experiment (shaded rectangles), with the corresponding uncertainties. The BFS results are 
indicated with arrows. 


always significantly greater than zero. The minimum value 
obtained is 54 pet for Ni^A^T^. In contrast, Cu(Al) 
decreases much more rapidly than Ti(A1) for increasing Al 
concentration and decreasing Ni and Cu concentration, even- 
tually becoming zero. This behavior indicates a switch in 
site preference exclusively favoring Ni sites (a Cu Ni = 0). 
This can also be seen in Figure 4(b), in the shape of a 
vaguely defined “boundary.” This transitional regime, 
loosely defined by alloys with a 1:1 Ni:Al ratio, is a region 
where P[Cu(Al)] (or the matrix element tf CuNi ) displays a 
sharp drop for decreasing Ni:Al ratios and where the corres- 
ponding Ti site preference for Al sites CPfTi(Al)], or a TlNi ) 
in Figure 4(a) decreases at a much slower rate and is mostly 
independent of Cu concentration. In this sense, mapping of 
the simulated results indicates much greater sensitivity of 
P[Cu(Al)] with concentration than that seen for P[Ti(Al)]. 
While the experimental results suggest a sudden reversal 
in site preference for Cu at a specific stoichiometry, the 
theoretical results indicate that this is a more gradual change 
beginning for alloys with 1:1 Ni:Al ratio (/.<?., the “surface” 
that could be built from the data displayed in Figure 4 for 
either F{Ti(Al)] or F[Cu(Al)] would resemble a “cascade,” 
smoothly evolving from one regime to another, as opposed 
to a sharp “step,” where there is no transitional region). 
We can also draw some general conclusions regarding the 


interaction between Ti and Cu atoms in the NiAl matrix. 
While a TlTl = 0 in every single alloy studied (indicating a 
strong repulsion between Ti atoms), a Cu cu is (in most cases) 
finite, allowing the possibility of clustering of a small frac- 
tion of Cu atoms, particularly for Ni-rich alloys (for example, 
a CuCu = 3.0 for Ni 5 oAl 43 Ti 3 Cu 4 , indicating that a Cu atom 
has a 3 pet probability of having another Cu atom as a NN). 
Moreover, the interaction matrix elements a TlC u and a CuTl are 
generally small for Ni-rich alloys, indicating the presence 
of few Ti and Cu NN pairs and, therefore, a higher likelihood 
of finding both elements in Al sites. These probabilities 
slowly increase as the concentration of Ni decreases, consis- 
tent with the fact that, for this range of concentrations, 
F[Cu(Al)] decreases much faster than P[Ti(AI)], thus 
favoring the location of Ti and Cu atoms in different sublat- 
tices. This results in a switch of site preference for Cu 
atoms from Al to Ni sites with a change in stoichiometry, 
as observed experimentally. 

Summarizing the results from these simulations, the strong 
ordering tendencies of NiAl alloys, coupled with the strong 
preference of Ti for Al sites, as well as the less prominent 
interaction between Cu and Ti atoms and the small -energy 
difference for Cu atoms in Ni vs Al sites, all contribute to 
make Cu atoms the ones most likely to fill in for any defi- 
ciency on either side of stoichiometry. This translates into 
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Table II. Coordination Matrices for the 5 NiAl-Ti-Cu 
Alloys Identified in Figure 1 and Displayed in Figure 3. 
Each Matrix Element a mn is Related to the Probability (in 
Percent) That an Atom n Has an Atom m as a NN, 
P[m(/i)]; for Example, the Probability That Ti has Ni as a 
NN, a r ,Nb atl Alloy of Composition NisoAl 41 Ti 3 Cu 6 is 
98.4 Pet, Indicating That Ti Resides Almost Exclusively on 
the Al Sublattice 


m\n 


Ni 

Al 

Ti 

Cu 

(a) Ni 5 oAl 41 Ti 3 Cu 6 

Ni 

4.3 

80.2 

6.0 

6.6 


Al 

97.7 

0 . 

0.0 

2.3 


Ti 

98.4 

0.0 

0.0 

1.6 


Cu 

80.3 

15.6 

0.8 

3.3 

(b) NisoAWTLC^ 

Ni 

2.7 

86.6 

6.0 

4.7 


Al 

98.5 

0.0 

0.0 

1.5 


Ti 

99.2 

0.0 

0.0 

0.8 


Cu 

77.0 

21.4 

0.8 

0.8 

(c) Ni 5 oAl 46 Ti 3 Cuj 

Ni 

1.2 

91.4 

6.0 

1.3 


Al 

99.4 

0.0 

0.0 

0.6 


Ti 

99.6 

0.0 

0.0 

0.4 


Cu 

68.8 

27.5 

1.2 

2.5 

(d) Ni49Al 47 Ti 3 Cui 

Ni 

0.0 

93.7 

6.1 

0.2 


Al 

97.8 

0.3 

0.0 

1.8 


Ti 

98.4 

0.8 

0.0 

0.08 


Cu 

10.0 

87.5 

2.5 

0.0 

(e) Ni 4 7 Al 47 Ti 3 Cu 3 

Ni 

0.0 

93.8 

6.2 

0.0 


Al 

93.8 

0.0 

0.0 

6.2 


Ti 

96.4 

0.0 

0.0 

3.6 


Cu 

0.0 

96.4 

3.6 

0.0 


a smooth transition from Al to Ni site occupancy as the 
change in composition becomes Ni poor. 

3. BFS analytical calculations 

Large-scale simulations provide, at best, a theoretical 
confirmation of the experimental results and, through the 
coordination matrices, some general insight on the possible 
underlying trends. Monte Carlo simulations provide useful 
information regarding the most likely final state of a given 
system. However, BFS-based detailed analytical calcula- 
tions can provide information on the structure of the low- 
lying energy states (i.e., most likely to appear) and the 
energetics of individual atoms or groups of atoms. A large 
set of reasonably large atomic configurations, each dis- 
playing a different distribution of atoms in a given computa- 
tional cell, is defined and its energy of formation computed 
using BFS. In this work, it proves to be sufficient to work 
on a 72-atom cell, as most relevant atomic distributions can 
be properly described. In what follows, we will use the 
notation Ni^Al^Ti^CUf/] to denote the concentration of the 
alloy in terms of the number of atoms in the 72-atom cell. 
Therefore, [i] = 36 corresponds to 50 at. pet Ni. 

In order to match the compositions studied either experi- 
mentally or with Monte Carlo simulations as closely as 
possible, we define 25 compositions that properly cover the 
entire range of Ni, Al, and Cu concentrations studied before. 
These states, denoted with gray solid circles in Figure 1, 
correspond to alloys Ni [A jAl [B jTi [2 jCu [C }> where the subindex 
indicates the number of atoms of each species in the 72- 
atom cell (A + B + C = 70). For each concentration, we 
build a catalogue of configurations that, if large enough, 
contains every possible arrangement of atoms that is likely 
to occur in the real alloy. We then compute the energy of 




Fig. 4 — Matrix elements (a) <z TlNi and ( b ) a Cu Ni> taken as a measure of 
P[Ti(Al)] and P(Cu(Ai)], respectively, as a function of Al and Cu concentra- 
tion, for different values of x Si . The projection of these curves onto the 
horizontal plane is shown with dashed lines, highlighting the fact that Cu(Al) 
becomes zero for high Al concentration and low Ni and Cu concentration. 

formation of each cell and plot the results in the form of 
an energy level spectrum, where, for each concentration, we 
can understand what atomic configurations are energetically 
favored and which ones are not. The number and type of 
configurations included in each catalogue are, in principle, 
arbitrary. In this work, a total of 300 configurations were 
used (for the entire range of concentrations studied). Exam- 
ples of configurations included in the catalogue are shown 
in Figure 5, where some selected configurations for one 
particular concentration, Ni^AlpojTi^Cu^, are shown. 
These cells can be described by labeling their differences 
with a basic B2 ordered NiAl cell. It is therefore convenient 
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Fig. 5 — (a) through (m) Some configurations for Nip^AlpojTipjCu^. Based on a NiAl 72-atom cell where the Ni sublattice is denoted with solid lines and 
the A1 sublattice with dashed lines, Ti additions are denoted with cyan (blue) disks if located in the Ni (Al) sublattice, and Cu additions are noted with 
magenta (red) disks if located in the Ni (Al) sublattice. Antistructure Ni atoms are denoted with large black disks. Otherwise, small black (gray) disks 
denote Ni (Al) atoms in their respective sites. The arrow indicates the first atom listed in the expressions displayed in Table III, where a general description 
of the chain of substitutional defects is provided. 
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(b) X(Ti)Ti(Ni) j Ni( Al) j > 2 Al(Ni) 1)2 ,i 


• Ni # Al O Ti X 

Fig. 6 — A Ni 2 AlTi cell after the introduction of an X atom in a Ti site 
(X(Ti)), the Ti atom moving to a Ni site (Ti(Ni)) so that the X and Ti 
atoms are located at NN distance, the Ni atom moving to an Al site (Ni(Al)) 
so that Ni is at NN distance with Ti and is a second neighbor of X, and 
the Al atom moving to a Ni site (Al(Ni)) so that it is NN of the Ni atom, 
second neighbor of the Ti atom and first neighbor of the X atom. 

to describe the relative location of each atom in the sample 
by means of a simple notation, which, while by no means 
complete, provides a first glimpse at the structure of each 
cell. In what follows, A and B represent the two simple 
cubic sublattices of the B2 compound and X represents a 
ternary alloying addition. The term X(A) denotes an atom 
X substituting for an atom A on the A sublattice, and 
X(A)A(B)„ denotes the same case but with the displaced A 
atom occupying a site in the B sublattice. The distance 
between atoms X and A is denoted by the subindex n. If 
n = 1 , atoms X and A are NN, while n = 2 indicates NNN. 

This notation can be easily generalized to describe more 
complex situations. In a chain of point defects, the first 
subindex denotes the distance (measured in NN or NNN 
distances) between the indexed defect and the one preceding 
it, and subsequent indices denote the separation of the 
labeled atom with previous atoms (from right to left) in the 
chain. For example, X(A)A(B) 1 B(A) 1)2 indicates an X atom 
occupying an A site and the A atom in a B site, so that X 
and A are NN (X(A)A(B)!). The B atom then moves to 
another A site, so that A and B are NN (the first subindex 
in B(A) lt2 ) and B and X are NNN (the second subindex in 
B(A) 1>2 ). Due to the fact that, in this application of BFS we 
deal with, at most, NNN interactions, we will not specify 
the distance between point defects if it exceeds the NNN 
distance. As an example, Figure 6 shows a variation of the 
perfect Ni 2 AlTi (Ll 2 ) Heusler cell where a chain of point 
defects is generated by the introduction of an atom X in a 
Ti site. The Ti atom in that site moves to a neighboring Ni 
site and, in turn, the displaced Ni atom moves to an Al site 
and, finally, the Al atom moves to an available Ni site. The 
X and Ti atoms are located at NN distance, the Ni and Ti 
atoms are located at NN distance and the Ni and X atoms 
are at NNN distance. Finally, The Al atom is a NN of the 
Ni atom, a NNN of the Ti atom, and a NN of the X atom 


(subindices 1, 2, and 1, respectively). The series of defects 
then reads X(Ti)Ti(Ni) t Ni( Al) ! 2 AI(Ni ) ! , 2 , i - In an expression 
dealing with several defects, the symbol X(A)_„ (n = 1 
and 2) means that atom X is at NN (n = 1) or NNN (n - 
2) distance from the first atom listed in the expression. 

We now return to the sample catalogue shown in Figure 
5 and write the corresponding series of substitutions using 
this notation in Table III. After minimizing the energy of 
formation of each cell with respect to the lattice parameter, 
we can order the configurations in an energy level spectrum, 
together with the resulting value of the lattice parameter, 
as shown in Table III. We also show the number and type 
of substitutional defects noting, not surprisingly, that the 
states with lower energy are those where Ti and Cu occupy 
Al sites. 

4. Ground state analysis 

Having shown the results for the ground state search for 
one specific concentration (Ni [36] Al [3 o]Ti r2 iCu [4] ), we now 
repeat this procedure for all the remaining compositions 
noted in Figure 1 and determine, in each case, the ground 
state. First, we concentrate on the structure of the ground 
states as a function of Cu concentration. The strong site 
preference of Ti for Al sites guides this choice, as it allows 
for a clear understanding of the behavior of Cu as it increases 
its role in the system. A series of Ni [38 _ t )Al [32+ ;c}Ti [2 ], (x = 
1, . . . , 5) alloys helps establish the behavior of Ti, when it 
is the only alloying addition. The behavior is dictated by 
the strong site preference of Ti for Al sites, regardless of 
the ratio between Ni and Al atoms. This preference for Al 
sites is absolute, in that Al antistructure atoms are created, 
if necessary, to accommodate all the Ti atoms in the Al 
sublattice (Ti(Al)), as shown in Figure 7. The Ti(Al) atoms, 
however, are somewhat sensitive to the presence of antisite 
defects, whether these are Ni atoms in the Al sublattice 
(which attract other Ti(Al) atoms at NNN distance) or Al 
atoms in the Ni sublattice (which attract Ti(Al) atoms at NN 
distance). The number of antistructure Al atoms is reflected 
in the rapidly increasing energy of formation for this set of 
alloys, with an average increase of 0.055 eV/atom per antisite 
defect. The actual atomic distribution in each alloy, as well 
as a simple schematic representation of these states high- 
lighting the behavior of the added Ti atoms, is shown in 
Figure 7. 

The Nipg-^Alfsj+^Ti^CUfn, (x = 1, .... 5) alloys are 
shown in Figure 8, following the same format and notation 
used in Figure 7. This set of alloys shows the first indication 
of the interaction between Ti and Cu. Once again, Ti(Al) 
substitutions dominate, leaving Cu second choice for avail- 
able Al sites. The Ti(Al) preference is strong enough as to 
induce the creation of Al antistructure atoms. Moreover, 
Ti(Al) seems to attract available Cu(Al) atoms along the 
[100] direction, thus inducing an ordered pattern where Ti 
and Cu atoms occupy alternating sites in the Al sublattice. 
However, the interaction between Ti(Al) and Cu atoms 
seems to be restricted to those cases where Cu goes to Al 
sites only. If Cu atoms are forced to go to Ni sites, they 
ignore the presence of Ti(AI) linking themselves only to Al 
antistructure atoms. These results highlight dominant (Ti(Al) 
inducing antistructure atoms) and secondary (Cu(Al) linked 
to Ti(Al)) features characterizing this group of alloys. So 
far, the main features that characterize the N Cll — 0 and N Cu 
= 1 cases are (1) absolute preference of Ti for Al sites; (2) 
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Table III. Description of the Cells Displayed in Figure 5 in Terms of the Different Substitutional Atoms, Relative to a 
Perfectly Ordered B2 NiAi Celi; the Number of Substitutional Atoms in Each Expression, the Difference in Energy 
of Formation (in eV/Atom) between the Cell and a Perfectly Ordered B2 NiAI Cell and the Value of the Lattice Parameter 

(in Angstroms) That Minimizes Such Difference Are Also Shown 


Cell 

Description 

Ti(Al) 

Ti(Ni) 

Cu(AI) 

Cu(Ni) 

Ni(Al) 

AH-AH 0 

a 

a 

Cu(Al) + Ti(Al)? + Cu(Al), + Ti(Al) 2 + Cu(AI) 2f2 + Cu(Al) f2 

2 

0 

4 

0 

0 

0.04274 

2.8448 

b 

Cu(Al) + Ti(Al) 2 + Cu(AI) 2 + Ti(AI) 2 + Cu(A1) 2 + Cu(Al) 

2 

0 

4 

0 

0 

0.04460 

2.8450 

c 

Cu(Al) 4- Cu(Al) + Ti(AI) + Cu(Al) + Ti(Al) 2 + Cu(Al) 

2 

0 

4 

0 

0 

0.04578 

2.8451 

d 

Cu(Al) + Ti(Al), + Cu(A1) 2 + Ti(Al) + Cu(AI) + Cu(Al) 

2 

0 

4 

0 

0 

0.04588 

2.8451 

e 

Cu(Ai) + Cu(A1) 2 + Ti(AI) 2 + Cu(A1) 2 + Ti(AI) 2f2 + Cu(Al) f2 

2 

0 

4 

0 

0 

0.04796 

2.8454 

f 

Cu(Al) + Ti(AI), + Cu(Al), + Cu(AI), + Cu(Al), + 
Ti(Ni), + Ni(AI),. a 

1 

I 

4 

0 

1 

0.12831 

2.8517 

g 

Cu(Al) + Ti(Al) 2 + Cu(Al) 2 + Cu(AI) 2 + Ni(AI) 2 . w + 
Ti(Ni)f tI j + Cu(Al) lif(2 

1 

1 

4 

0 

1 

0.15883 

2.8546 

h 

Ti(Ni) + Ti(Ni) 2 + Cu(AI) M + CuCAl)^,, + Ni(Al) p . 2 , M + 
Cu(Al) 2 ,fxi + Cu(A1) 2 j + Ni(Al) 2- f, 2 , 2 ,f,i.i 

0 

2 

4 

0 

2 

0.18520 

2.8565 

i 

Cu(Ni) + Nfi(Al) | + Ti(Ni), + Ni(Al), + Cu(Ni), ff2 + 
Ni(Al)i j + Ti(Al) 2 , + Cu(Ni), + Ni(AI), + ’ 
Cu(Ni).:. xu + Ni(Ai)! 

1 

1 

0 

4 

5 

0.22614 

2.8609 

j 

Cu(Ni) + Cu(Al), + Ti(Ni), 4- Cu(Ni) , + Ni(Al), + 
Ni(Al) f , j + Cu(Al) 22 i + Ti(Ni)i + Ni(Al), , + Ni(AI) 

0 

2 

2 

2 

4 

0.27844 

2.8651 

k 

Cu(Ni) 4- Ni(AI), + Ni(Al) 2I + Ni(AI), fI 4- Ti(Ni), j + 
Cu(Ni) f , 4- Cu(Ni) 2 4- Ni(Al), + Ni(Al), , + Ti(Ni), 
+ Cu(Ni) + Ni(Al)! 

0 

2 

0 

4 

6 

0.28864 

2.8661 

I 

Ti(Ni) + Cu(Al)j + Cu(Al) 2 + Cu(A1) 2 + Ni(Al) 2m + 
Cu(A1) 2 + Ni(Al) 2 ,f, 2 4- Ti(Ni)! 

0 

2 

4 

0 

2 

0.31041 

2.8682 

m 

Cu(Ni) + Ni(Al), + Cu(Ni) 1 4- Ni(Al) ifl + Ti(Ni)j 4 Ni(AI), 
+ Cu(Ni), + Ni(Al), Ni(Al) f , f , + Ti(Ni) 4 Cu(Ni) 2 4 
Ni(Al) 1(1(f>ftl 4 

0 

2 

0 

4 

6 

0.32964 

2.8701 


first choice of Ti for A1 sites, overriding Cu’s preference 
for A1 sites; (3) creation of antistructure A1 atoms when the 
number of Ti atoms exceeds the number of available A1 
sites (i.e. y for N Al >34); (4) clustering of Ti(Al) around 
Al(Ni) antisite defects; (5) clustering of Cu(Al) around 
Ti(Al) (at NNN distance); and (6) decoupling between Ti(Al) 
and Cu atoms when Cu goes to available Ni sites. 

The iV Cu = 2 (Ni l 3 8 ^^]Alpo+.xiTi[ 2 jCuf 2 ], (x — 1, . . . , 5)), 
shown in Figure 9, provides information on the interaction 
between Cu atoms and the enhanced competition between 
Ti and Cu additions with increasing Cu concentration. A 
common feature of the alloys in this set continues to be the 
absolute preference of Ti for A1 sites and the ensuing creation 
of antistructure A1 atoms when A1 sites become unavailable. 
The Ti(Al) clusters form around the Ni antistructure atom. 
If A1 sites are available, Cu(Al) atoms couple to Ti(Al) 
atoms, but as more Cu moves to the Ni sublattice, the 
coupling with Ti(Al) is replaced by either complete 
decoupling (as in Nip^Al^Ti^jCupj) or attraction with 
other Cu atoms (as in Ni^AlpsjTi^jCu^). The /V Cu = 3 
group (Ni [38 _. c] Al [2 9 + ^Tir 2 l Cu [3 |, for x = 1, .... 5), shown 
in Figure 10, highlights the behavior of Cu, as its percentage 
now exceeds that of Ti. The only additional feature observed 
in this group, a consequence of the increased Cu content, 
consists of the apparent clustering of Cu atoms (either in 
the Ni or A1 sublattice), decoupled from any Ti(Al) atoms. 
The diagrams in Figure 10 succinctly show the competition 
between Ti-Cu ordering in the A1 sublattice and the cluster- 
ing of Cu atoms in either sublattice. The N Cll — 4 group 
(Ni| 38 - rj Alps +^]Ti f2] Cu| 4 j, for* = 1, , 5), shown in Figure 
1 1, displays trends already apparent from the previous cases. 
Although no new features are observed, the larger number 


of Ti and Cu atoms helps to understand the main characteris- 
tic of the Ni:Al 1:1 alloys, where most additions are in 
solution and the clustering tendencies clearly observed in 
other systems are much less pronounced (i.e., coupling of 
Ti(Al) and Cu(Al) atoms or clustering of Cu atoms in either 
sublattice). In all cases, as expected, the minimum energy 
occurs for Ni [rt] Al [/? |Ti [2] CU[ c j alloys, where b + c + 2 = a, 
and where no antistructure atoms exist. In terms of changes 
in the lattice parameter, it is observed that in all cases the 
increase in A1 concentration is the only source of such 
changes, due to the small difference in size between Ni and 
Cu atoms. 

It is interesting to group the results, individually shown 
in Figures 7 through 1 1 , in one single figure as a function 
of Ni, Al, and Cu concentration, as shown in Figure 12. 
Visual examination provides an indication of the characteris- 
tic features of each region in the range of concentrations 
studied: the Ti-Cu ordering in the Al sublattice that character- 
izes alloys where N m > N A] transitions to an intermediate 
regime where no specific pattern dominates (V Ni ~ V A1 ), to 
be replaced later by another ordering pattern in the Ni sublat- 
tice (for Al-rich alloys) governed by Cu-Cu interactions. 

IV. DISCUSSION: LOCAL ENVIRONMENT 
APPROACH 

The picture that emerges from Figure 12, while retaining 
the essential features already observed in Monte Carlo 
simulations, describes the results of the interaction between 
different alloying additions, Ti and Cu, as a function of 
composition. To gain a proper understanding of this interac- 
tion, it is useful to perform an atom-by-atom energy analysis 
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a. Ni [37 jAl [ 33jTi [2 ] 

- 0.57937 / 2.8536 
Ni(Al)+Ti(Al)2+Ti(Al)3 )2 
Excess Ni(Al), two direct 
Ti(Al) substitutions 


b * NipflAl[34jTI(2] 

- 0 . 58580 / 2.8580 
Ti(Al)+Ti(Al) f 
Two direct Ti(Al) in 
solution 



e. Nip 5 ]AI{35]Tl[2] 

- 0.52298 / 2.8690 
Ti(Al)+Ti(Al)Al(Ni) u 
One direct, one ‘forced’ 
Ti(Al), linked by Al(Ni) 



Pm 


d. Ni[34]Al[36]T^[ 2 ] 

- 0.46710 / 2.8793 
Ti(Al)Al(Ni)] +Ti(Al)Al(Ni)j 2 
Two ‘forced’ Ti(Al) substitutions, 
linked by both Al(Ni) atoms 


ra 


Ti 

— | 

Al 

T 

i® 

I 

M 


[Al 

] : 

Ti 

M 

Al 



e. Ni [33 ]Al [ 37]Tl [ 2] 

- 0 . 41384 / 2.8900 
A1(N i)+Ti(Al)Al(Ni) i 2 
+Ti(Al) 1 Al(Ni) 1>2 j 

Excess Al(Ni), two forced Ti(Al) 
atoms linked by Al(Ni) atoms. 




Fig. 7 — ( a ) through ( e ) Ground state configurations including the energy of formation (in eV/atom), lattice parameter (in Angstroms) and structural 
information for Nipg-^Al^+jjTi^, for x = 1, . . . , 5. Squares denote atoms in the AI sublattice, while circles denote atoms in the Ni sublattice. The arrows 
in (d) and (e) indicate whether the atom is in a plane above (up) or below (down) the plane of the page. On the right, a simpler, schematic diagram 
summarizes the principal features of the atomic distributions shown on the left, concentrating only on the defect structure. Solid vertical or horizontal Tines 
indicate bonds in an Al plane, while solid diagonal lines (c, d, e) indicate bonds between an atom in the Al sublattice and its NN in the Ni sublattice. 
Dashed vertical or horizontal lines indicate NNN bonds between atoms in the Ni sublattice. Al* denotes an antistructure Al atom, Al(Ni), and Ni* denotes 
an antistructure Ni atom, Ni(Al). 


in order to identify the trends and reasons for the observed 
behavior. A simple approach, based on the BFS energy con- 
tributions for a given defect and its surroundings, helps 
explain the observed behavior in Ni-Al-Ti-Cu alloys, as well 
as the magnitude of the energy gaps that characterize the 
various defect structures observed in the energy spectra of 
the alloys studied. Based on the energetics of small groups 


of atoms surrounding specific atoms in the cell, the idea 
consists of analyzing the resulting site substitution behavior 
not just in terms of individual interactions between individ- 
ual atoms, but through the effect that these bonds might 
have in the immediate vicinity of such atoms. 

In order to implement this approach, we define atomic 
“local environments” (LEs) fl8] consisting of a given central 
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a. Ni[ 3 7 ]A.l[ 32 ]Ti[ 2 ] Cu [i] 
- 0 . 57330 / 2.8502 

Ni(Al)+Ti(Al) 2 +Cu(Al) 2 +Ti(Al) 2 , fi2 
Excess Ni(Al), two direct Ti(Al) and 
one direct Cu(Al) substitutions. 


b - Ni[ 36 ]Alp 3 ]Ti [ 2 jCu [1 ] 

- 0.58227 / 2.8545 

Ti(Al)+Cu(Al) 2 +Ti(Al) 2 

Two direct Ti(Al) and one Cu(Al) substitutions, 

linking both Ti(Al) atoms. 



c - Ni [35] A, [34] Ti [2] Cu [l] 

- 0 . 57239 / 2.8608 
Ti(Al)+Ti(Al)f+Cu(Ni) f f 
Two direct Ti(Al), one direct Cu(Ni) 
substitution in solution. 




d * Ni [34] A, [35] Ti [2] Cu [l] 

- 0 . 51152 / 2.8714 

Ti( Al)+Ti ( Al) Al(Ni) 1 +Cu(Ni) 2 

One direct, one ‘forced’ Ti(Al) substitutions, 

linked by both Al(Ni) atoms, and one direct 

Cu(Ni) substitution linked to Al(Ni), but not to 

Ti(Al). 




e - Ni^Al^jTipjCUfi] 

- 0 . 45967 / 2.8834 

Ti(Al) Al(Ni) i +Ti(Al) Al(Ni) 12 +Cu(Ni) 2 
Two forced Ti(Al) atoms linked by Al(Ni) atoms. 
One direct Cu(Ni) linked to Al(Ni) but not to 
Ti(Al). 


Fig. 8 — (a) through ( e ) Ground state configurations for Ni f 3 8 _ Jcl AI r3 , ^ xl Tij- 21 Cu r , for x - 1, . . . , 5. 



atom and its eight NN and six NNN (in a bcc lattice), under 
the assumption that it is this group of atoms that will be 
most affected by the presence of a central substitutional 
defect atom and that any change in the energy of formation 
of the cell will arise mostly from changes in energy within 


this environment. Considering the fact that the BFS equa- 
tions deal with up to second neighbors of a given atom, and 
that no individual relaxations are allowed, it is then sufficient 
to examine the energetics of this limited group of atoms in 
order to understand the full effect of a given point defect 
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a. Ni[37]Ai[3i]Tip]Cu[2] 

- 0.56688 / 2.8468 

Ni(Al)+Ti(Al) 2 +Ti(Al) f>2 +Cu(A 1) 2 + Cu(Al) 
Excess Ni(Al), two direct Ti(Al) and two direct 
direct Cu(Al) substitutions. 



b - N5 [36} A1 [32] Ti [2] Cu [2] 

- 0 . 57529 / 2.8512 

Ti(Al)+Cu(Al) 2 +Ti(Al)2 *fCu(Al)2 f 2 
Two direct Ti(Al) and Cu(Al) 
substitutions. 



c - Ni [35] AI [33] Ti [2] Cu [2] 

- 0 . 56617 / 2.8574 
Cu(Ni)+Cu(Al) 1 +Ti(Al) ? +Ti(Al) 
Two direct Ti(Al), one direct Cu(Ni) 
linked to Cu(Al). 


© 



d.Ni[34]Al t 34 ] Ti[ 2 ]CU[2] 
- 0 . 55662 / 2.8636 
Ti(Al)+Ti(Al)+Cu(Ni) + Cu(Ni) 

Two direct Ti(Al) substitutions and two 
direct Cu(Ni) substitutions, in solution. 



- 0 . 49871 / 2.8740 

Ti(Al)Al(Ni) 1 +Ti(Al) 1 +Cu(Ni)+Cu(Ni) ? 

One direct and one forced Ti(Al) atom linked by 
Al(Ni). One direct Cu(Ni) linked to Al(Ni) but not 
to Ti(Al). 



Fig. 9 — (a) through ( e ) Ground state configurations for Ni I38 _ r] Al r 3 0 ^ r| Ti (2I Cu f21 , for.r = I, .... 5. 


(/.£., substitution of one atom by another, creation of a single 
vacancy, etc.). As an example, the LE of an individual X(B) 
defect, denoted by (X(B)), is shown in Figure 13. The (XQB)) 
consists of an X atom in a site in the B sublattice of an A-B 
alloy, surrounded by eight NN of type A and six NNN of 
type B. The formation of an antisite defect, X(A)A(B), could 


then be seen as the superposition of the LE surrounding two 
individual point defects: one centered around the X atom, 
(X(A)), and the other centered around the displaced A atom, 
(A(B)). These could be “noninteracting.” ([X(A)A(B)] / ) (if 
the two point defects are separated by more than one lattice 
parameter distance), or “overlapping.” ([X(A)A(B)]i), LEs, 
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a - Ni [37] AI [30] Ti [2] Cu [3] 

- 0.55866 / 2.8437 

Ni(Al)+Ti(Al) 2 +Cu(Al) 2 +Ti(Al) 2 +Cu(Al) 2>fi2 
+Cu(Al)f 2 

Excess Ni(Al), two Ti(Al) and three direct 
Cu(Al) substitutions 


b - Ni [36] A1 [31] Ti [2] Cu [3] 

- 0.56789 / 2.8480 

Ti(Al)+Cu(Al) 2 +Ti(Al) 2 +Cu(Al) 2f2 + Cu(A 1) 2 
Two direct Ti(Al) and three direct Cu(Al) 
substitutions 


c* Ni [35] A1 [32] Ti [2] Cu [3] 

- 0.55847 / 2.8541 

Ti(Al)+ Cu(A1) 2 +Cu(A 1) +Cu(Ni) 1 +Ti(Al) 
Two direct Ti(Al), one direct Cu(Ni) linked to 
Cu(Al) 


d. Ni[34]Al[33]Ti[2]Cu[ 3 ] 

- 0.55085 / 2.8602 

Ti( Al)+Ti (Al)+Cu(Ni) + CuCAlA+CuCNi)! 

Two direct Ti(Al) substitutions, one Cu(Al) and 
two Cu(Ni) substitutions, with clustering of Cu 
atoms. 


e * Ni [33] A1 [34] Ti [2] Cu [3] 

- 0 . 54248 / 2.8663 

Ti(Al) 1 +Ti(Al)+Cu(Ni)+Cu(Ni) 2 +Cu(Ni) 2 
Two direct Ti(Al) atoms in solution, a cluster of 
three direct Cu(Ni) atoms, not coupled to any 
Ti(Al) atom 



Fig. 10 — (a) through {e) Ground state configurations for Ni f 3 8 _ j] Al (2 9 ^ c1 Ti r2| Cup,, for jt = 1, . . . , 5. 


where the two center atoms (X and A) are NN, thus sharing To make a fair comparison of the different single and 

a large number of NN and NNN atoms. The latter case is extended local environments, it is necessary to embed each 
shown in Figure 13(b). LE in a B2 cell so that all bonds affected by the presence 
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a. Ni [37] Al [ 29]Ti [ 2jCu [ 4] 

-0.55125/2.8406 

Ni(Al)+Ti(Al)+Cu(Al) 2 +Ti(Al) 2 +Cu(A 1) 2 + 
Cu(Al)f 2 

Excess Ni(Al), two Ti(Al) and four direct 
Cu(Al) substitutions. 


k Ni[ 36 ]Alpo]Ti[ 2 ]Cu [4] 

-0.56030/2.8450 

Ti(Al)+Cu(Al) 2 +Ti(Al) 2 +Cu(Al) 2 f>2 + Cu(A 1). 2 
Two direct Ti(Al) and four direct Cu(Al) substitu 
tions. 


Ni [35] A, [31] Ti [2] Cu [4] 

-0.55338 / 2.8507 

Ti(Al)+ Cu(Al) 2 +Ti(Al)+Cu(Al) 2 f2 
+Cu(Ni) 1 +Cu(Al) ] 

Two direct Ti(Al), one direct Cu(Ni) linked 
to Cu(Al). 


d* Ni[ 34 ]Al t 32 ]Ti [ 2 ]Cu [4 j 
-0.54786/2.8564 

Ti(Al)+Cu(Al) 2 +Ti(Al) 2 +Cu(Al) 22 +Cu(Ni) 1 +Cu 

(Ni) 2 

Two direct Ti(Al) substitutions, two Cu(Al) and 
two Cu(Ni) substitutions, with clustering of Cu 
atoms. 


e * Ni [33] A1 [33] Ti [2] Cu [4] 

-0.53742/2.8628 

Ti(Al)+Cu(Al)2+Cu(Ni)l+Cu(Ni) 2(] +Ti(Al)+Cu(Ni) 
Two direct Ti(Al) atoms in solution, one Cu(Al) 
and three direct Cu(Ni) atoms, not coupled to any 
Ti(Al) atom. 



Fig. 11 — (a) through ( e ) Ground state configurations for Ni^g^Al^g+jrjTi^Cu^, for x = 1 5. 


of the defect are accounted for. To do so, we locate the LE 
at the center of a 72-atom equilibrium B2 NiAl cell and 
refer all energies to the pure AB version of that cell. Because 


all LEs, as well as the reference B2 cell, are evaluated at 
the same lattice parameter, the energy difference between 
the cell with the defect and the reference cell, Ae, indicates 
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• •a Oa(B) ©X • •B 

Fig. 13 — {a) (X(B)) local environment: X atom (light gray) in a B site, 
surrounded by eight A NN (large black disks) and six B NNN (large gray 
disks). ( b ) <[X(A)A(B)]i), with the substitutional X atom in an A site and 
the displaced A atom (circle) in a B site. The two “defect” atoms (the 
substitutional X and antistructure A) are NN and are surrounded by an 
environment of 20 atoms, accounting for all the NN and NNN, some of 
which are shared between them. A or B atoms at distances greater than 
one lattice parameter from either one of the defect atoms are indicated with 
small black and gray disks, respectively. 


the energy cost (in terms of chemical energy differences) of 
performing the specific substitutions that characterize the 
defect (as it contains only the contributions from the atoms 
in the LE, and not the rest of the cell). 

Table IV lists the energy of formation (relative to an ideal 
B2 NiAl cell) of different local environments: (X(B)) and 
(X(A)) (for X = Ti or Cu and A and B = Ni or Al), (Ni( Al)> 
and (Al(Ni)}. The corresponding values for the extended 
local environments, including two substitutions, are also 
listed. It is reasonable to expect that the site preference is 
mainly determined by the relative energy of formation 
Ae{LE} of these local environments. For example, if 
A*{(X(B)>} < A^{([X(A)A(B)] / >} or A*{<X(B)>} < 
A^{([X(A)A(B)]i)}, then X chooses a site in the B sublattice, 
as the formation of a single (X(B)> entails a lower energy 
cost than that required by the combination of a substitutional 
defect and the creation of an antistructure A atom, regardless 
of the relative position of these defects. For consistency, this 
analysis has to be performed in unrelaxed atomic positions, 


Table IV. BFS Energy of Formation of Single and Double 
Local Environments (in eV/Atom), Relative to a 72-atom 
B2 Cell; for the Double-Centered LE, the Center Atoms 
Are Located at NN ([] x ), NNN ([] 2 ), or Greater Than 
NNN Distance ([]/); the Subindex H ([]^) Indicates That the 
Center Atoms (Ti and Cu) Locate Themselves Following 
a Heusler (L2j) Pattern in the Al Sublattice (i.e., Sharing 
an Al Atom as a NNN) 


LE 

Ae T 

(Ni(Ni)> 

0.00000 

(Ni(Al)) 

0.01099 

(AlfNi)) 

0.09434 

(Ti(Al» 

0.00772 

(Ti(Ni)) 

0.16333 

(Cu(Ni)) 

0.01651 

(Cu(AI)) 

0.00889 

<[Ti(Al) + AI(Ni)]r) 

0.10205 

<[Ti(Al) + Al(Ni)],} 

0.08739 

<[Ti(Al) + Cu(AI)] 2 > 

0.01497 

<[Ti(Al) + Cu(Al)]„) 

0.01622 

<[Ti(Al) + Cu(Al)],) 

0.01661 

([Ni(Al) + Ti(Al)jf) 

0.01870 

<[Ni(AI) + Ti(Al)],> 

0.01703 


LE &e T 

<[Cu(Al) + Ni(Al)] f > 0.01988 
([Cu(Al) 4- Ni(Al)] 2 ) 0.02151 
<[Al(Ni) + Cu(Al)] f ) 0.10323 
([AJ(Ni) + Cu(Al)],> 0.09369 
([Ti(Al) + Ti(Al)] f ) 0.01544 
<[Ti(Al) 4- Ti(Al)] 2 > 0.01647 

([Ti(Al) + Ti(Al)] H > 0.01593 
([Cu(Al) + Cu(Al)] f ) 0.01778 
([Cu(AI) 4- Cu(Ni)] f ) 0.02540 
<[Cu(Al) + Cu(Ni)]!> 0.02377 
<[Ti(Ni) + NKADh) 0.15060 
([Ti(Ni) 4- Ni(Al)] f ) 0.17431 
([Ni(Al) 4- Cu(Ni))]) 0.02598 
([Ni(Al) 4- Cu(Ni)] f ) 0.02750 


under the assumption that relaxation effects, while influen- 
tial, are not ultimately responsible for the observed site 
preference behavior. 

A. Local Environment Analysis of Atomic Coupling 

In this approach, it is as important to isolate individual 
patterns and establish a hierarchy of effects by means of 
their relevance as measured by the typical energies involved, 
as it is to first gain understanding of the behavior of individ- 
ual additions before examining the effect of their interac- 
tions. We start with the LE study of ternary additions (Ti or 
Cu) to NiAl, followed by the quaternary case (Ti and Cu 
interactions in NiAl). In the ternary cases, where the goal 
is to identify the reasons that explain the different site substi- 
tution options, it is not necessary to consider the relative 
locations of the substitutional and antistructure atoms. While 
impor-tant, the gain or loss of energy due to different relative 
positions is not comparable in magnitude with the corre- 
sponding gain or losses due to substitutions themselves. 

1 . 77 site preference in NiAl 

The relevant results listed in Table IV are shown in 
Figures 14(a) and (b), which schematically represent the 
energy cost (in terms of the energy of formation of unre- 
laxed LEs) in creating substitutional defects and anti struc- 
ture atoms for Ni-rich (Figure 14(a)) and Al-rich alloys 
(Figure 14(b)). 

For Ni-rich alloys, which contain available Al sites, the 
competition between Ti(Al) and Ti(Ni)Ni(AI) substitutions 
is explained by the high energy cost of the Ti(Ni) substitu- 
tion ((Ti(Ni))), which, compounded with the small energy 
cost of creating a Ni antistructure atom ((Ni(Al))), 
results in an almost insurmountable energy gap greatly 
favoring Ti(Al) substitutions: Ae{ (Ti(Ni)Ni(Al))/} - 
Ae{(Ti(Al))} = 0.16660 eV/atom. In Al-rich alloys (Ni 
available sites), it is the high energy cost of creating an 
antistructure Al atom that dominates the site preference 
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<Ti(Ni)Ni(Al)> 
Ti in Ni sites 


<Ti(Ni)> 
Ti in Ni sites 


| Ni(Al) 


Ti(Ni) 


<Ti(Al)> 
Ti in A1 sites 


x Ni > x M 

(a) 


<Ti(Al)Al(Ni)> 
Ti in A1 sites 


Al(Ni) 


f Ti(Al) 

x Ni < X A1 

(b) 


Ti(Ni) 


▼ 


<Cu(Al)Al(Ni)> 
Cu in Al sites 


<Cu(Al)> 
Cu in Al sites 


<Cu(Ni)Ni(Al)> 

Cu in Ni sites 

“■'""J Ni(Al) 

"" i Cu(Ni) 


Al(Ni) 


<Cu(Ni)> 

Cu in Ni sites 

r — 

Cu(Al) 


iCu(Ni) 


x Ni > X A1 x Ni < X A1 

(c) (d) 

Fig. 14 — Energy level diagram for the energies of formation of different local environments needed to describe the site preference behavior of (a) Ti in 
Ni-rich, ( b ) Ti in Al-rich, (c) Cu in Ni-rich, and (d) Cu in Al-rich NiAl alloys. The left diagrams indicate the possible substitutions in Ni-rich alloys, 
indicating a clear preference for Al sites for Ti and Cu. The diagrams on the right correspond to Al-rich alloys, still showing a clear (but less marked) Ti 
preference for Al sites and Cu preference for Ni sites. The energies are referenced to a pure B2 NiAl alloy. 


behavior, as seen in Figure 14(b). However, in spite of 
being large (A^{(Al(Ni))} = 0.09434 eV/ atom), the 
combination of a Ti(Al) and a Al(Ni) substitution still does 
not match the high energy cost of locating a Ti atom 
in a Ni site. The energy gap between Ti(Ni) and 
Ti(Al)Al(Ni) defects in Al-rich alloys, Ae{(Ti(Ni)}} - 
Ae { <Ti( Al) Al(Ni)}/} - 0.06127 eV/atom, is not as large 
as the one between Ti(Ni)Ni(Al) and Ti(Al) defects in Ni- 
rich alloys, indicating that, while the preference of Ti for 
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Al sites is shared by Ni-rich and Al-rich alloys, it is weaker 
for the latter, consistent with the results of the Monte Carlo 
simulations where there was an increasing proportion of 
Ti atoms finding their way to the Ni sublattice for Al- 
rich alloys. 

2. Cu site preference in NiAl 

The analysis of Cu additions to NiAl is shown in Figures 
14(c) and (d). For Ni-rich alloys (Figure 14(c)), the energy 
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cost of locating a Cu atom in a Ni site, added to a comparable 
energy cost for creating a Ni antisite defect, results in a net 
preference for AI sites, although not of the same magnitude 
as that found for Ti in a similar situation. The proximity 
between the energy levels corresponding to Cu(Al) and 
Cu(Ni)Ni(Al)/ can be understood as more options become 
available for Cu, which can be easily excited to reside in a 
Ni site with low energy cost. For Al-rich alloys, in spite of 
the fact of the direct preference of Cu for Al sites 
(A^{(Cu(Al))} = 0.00889 eV/atom vs Ae{(Cu(Ni)}} = 
0.0165 1 eV/atom), the high energy cost of creating an antisite 
Al defect erases the advantage of direct Cu substitutions for 
Al sites, making Ni sites a favorable choice. 


B. 77 and Cu Additions and Interaction between Point 
Defects 

Before continuing with a detailed analysis of the quater- 
nary case, some conclusions can be extracted from the ter- 
nary cases, as described in Figure 14. It is interesting to 
compare the Cu energy level diagrams with those obtained 
for Ti in spite of the fact that both elements display a “direct” 
preference for Al sites regardless of composition (i.e. 7 either 
Ti or Cu is energetically more stable in an Al site than in 
a Ni site). Comparing Figures 14(a) and (b) with Figures 
14(c) and (d) and the corresponding LE formation energies 
in Table IV, the energy difference between Cu(Al) and 
Cu(Ni)Ni(Al) substitutions is only 0.01861 eV/atom, an 
order of magnitude smaller than that found for Ti. This fact, 
coupled with the smaller energy cost in locating Ti atoms 
in Al sites in Ni-rich alloys than the much larger one required 
for substitutions in Ni sites, justifies the observed prevalence 
of Ti in using the available sites in the Al sublattice. Con- 
versely, for Al-rich alloys, the resulting energy difference 
between Cu(Al)Al(Ni) and Cu(Ni) substitutions (0.08672 
eV/atom) explains the ease with which Cu favors sites in 
the Ni sublattice instead of competing with Ti for Al sites. 
It could then be assumed that if this behavior continues in 
the case of quaternary systems, then Ti will always choose 
Al sites while Cu will move from Al sites to Ni sites 
depending on the Ni:AI ratio. This assumption, based solely 
on extrapolating features characterizing the ternary systems, 
ignores effects that could arise from the interaction between 
the two alloying additions. The atomic distributions shown 
in Figure 12 clearly indicate that the ternary site substitution 
behavior alone does not completely explain all the salient 
characteristics of each ground state, as several other features, 
besides the specified site preference for Ti and Cu, can be 
identified (Figure 12). 

It is necessary to concentrate on the quaternary system 
and the influence of the interaction between the different 
elements in the resulting structures. These results indicate 
that the characteristic energies of the single element site 
preference behavior (Ti or Cu in NiAl) are of the order of 
0.1 eV/atom. There are clearly secondary effects, at a lower 
energy scale, that also intervene in determining the final 
atomic distribution. Such effects arise from the coupling 
between the different alloying additions and should not be 
ignored if a proper description of the ground state is desired. 
Several of these effects have already been identified, 
although not explained, by inspecting the structure of the 


ground states shown in Figures 7 through 11. In what fol- 
lows, we describe and quantify the nature of these effects, 
based on the analysis of noninteracting and overlapping LEs. 

1 . 77 and Cu interaction with antisite defects 

One common feature of all the ground states shown in 
Fig. 12 consists of the preference of Ti for Al sites. It is 
interesting to discuss how Ti(Al) atoms interact with anti- 
structure atoms, when present, or Cu alloying additions. 

The top row in Figure 12 contains alloys for which it has 
been established that Ti and Cu additions occupy sites in 
the Al sublattice. Moreover, for x Ni > 50 at. pet, there are 
also Ni antistructure atoms present. With these three types 
of atoms in the Al sublattice, we now investigate the possible 
coupling schemes that might appear in such systems. By 
coupling, we mean the interaction between different types 
of atoms that results in the formation of NN or, when appro- 
priate, NNN bonds. 

Among the three possible couplings, [Ti(Al) 4- Ni(Al)], 
[Cu(Al) 4* Ti(Al)], and [Cu(Al) 4- Ni(Al)], only the first 
two involve a gain in energy when a NNN bond is formed 
(/.£., A^{[] 2 } < A<? {[]/}). The creation ofa[Ti(Ni) 4- Ni(Al)] 2 
bond introduces a gain of 0.00167 eV/atom over [Ti(Al) 4- 
Ni(Al)]/. This explains the observed Ti(Al) Ni(Al) cou- 
pling that characterizes the ground states of the Ni-rich qua- 
ternary alloys in the top row in Figure 12. A comparable 
gain is realized if a [Cu(Al) 4- Ti(Al)] 2 bond is formed. 
No such gain exists for [Cu(Al) + Ni(Al)] 2 bonds, thus 
establishing a coupling scheme Cu(AI) «-* Ti(Al) <-► Ni(Al). 
Table IV also shows that the local environment energy of 
Ti(Al) 4* Cu(Al) atoms varies little depending on the relative 
location of the two substitutional atoms (0.01497 eV/atom 
when they are located at NNN distance, 0.01661 eV/atom 
when they are separated by greater distances, and 0.01622 
eV/atom when they are located at the opposite corners of a 
cube in the Al sublattice). The proximity of these three 
energy levels explains the ordering pattern observed, for 
example, in the ground state for the Ni^Al^Ti^Cu^] 
alloy. 

The coupling scheme (Cu(Al) <-► Ti(Al) ++ Ni(AI)) results 
in the characteristic feature of other Ni-rich alloys, as shown 
in Figure 15. Some additional patterns with respect to that 
coupling hierarchy are also apparent. With the coupling 
Ti(Al) <-► Ni(Al) firmly established, we now consider its 
consequences when several Ti(Al) atoms are present. Of the 
three cases listed in Table IV for the relative location of two 
substitutional Ti(Al) atoms ([Ti(Al) + Ti(Al)]/, [Ti(Al) 4* 
Ti(Al)] 2 , and [Ti(Al) 4- Ti(Al)] H ), the NNN bond is the less 
energetically favored, followed by the “Heusler” arrange- 
ment and, finally, with the lowest energy cost, the two atoms 
far apart. As mentioned previously, Ti(Al) ++ Ni(Al) coupling 
will rule out this last option, thus favoring sites at third- 
neighbor distance (Heusler-like sites), as seen in Figure 15, 
where the number of favorable Ti-Al bonds is maximized. 
This effect can be understood as the role of Ni(Al) promoting 
the attraction of Ti atoms into an ordered pattern. While the 
link Cu(Al) <-> Ti(Al) can be interpreted as responsible for 
“closing” the square formed by Ni(Al), Ti(Al), Cu(Al), and 
Ti(Al) (Figure 15), the “repulsion” between Cu(Al) and 
Ni(Al) results in a competition between the number of 
Ti(Al)-Ni(Al) and Cu(Al)-Ti(Al) bonds, which explains the 
loss of one Ti(Al)-Ni(Al) bond at the expense of extra 
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Ni* — Ti + | | + Ti Cu 

Ti ' ' 

Cu coupling to Ti 
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Ti Ni* Ti 

1 1 
Cu Cu 

Ni* — Ti + Ti Cu 
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Cu 

1 
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| 

Ti Cu 
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Cu Ti Cu 

Ni* — Ti + Ti Cu 


Ni* 
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1 

Cu Ti Cu 

| | 

Abundance of Cu favors Ti/Cu coupling 
over Ti coupling to Ni 


Cu Ti Cu 



Fig. 15 — Coupling series Cu Ti «-*■ Ni, as manifested in Ni-rich alloys with increasing Cu content. The last column describes the individual elements, in 
order of importance, leading to the final state (center column). 


Cu(Al)-Ti(Al) bonds, as observed in Figures 15(d) and (e). 
However, the alternating pattern of Cu-Ti-Cu-Ti atoms per- 
sists. This fact can be taken as a hint that in alloys with 
higher Ti concentration than the ones studied in this work, 
Cu additions might partition to the L2 { (Heusler) phase. 

The coupling of Ti(Al) atoms with antisite defects persists 
in Al-rich alloys. For these cases, the Ti(Al)-Al(Ni) coupling 
introduces comparable energy gains (0.01466 eV/atom) with 
respect to the noninteracting case. Every alloy with A1 anti- 
structure atoms shows Ti(Al)-Al(Ni) coupling, regardless of 
the Cu contents. Moreover, Ti(Al)-Al(Ni) coupling overrides 
the tendency of Ti(Al) atoms to locate themselves in Heusler 
sites. Overall, this coupling seems to be the leading effect 
relating the defects in each sublattice in these alloys: Ti(Al) 
substitutional atoms in the A1 sublattice and Cu atoms in 
the Ni sublattice both bond to the available antistructure A1 
atoms (Figures 7(c), (e), 8(d) and (e), and 9(e)). 


2. Ti and Cu interactions 

As the number of available A1 sites decreases (W Ni < 36 
in Figure 12), the leading role of Ti(Al) substitutions forces 
Cu atoms to occupy the remaining A1 sites as well as the 
increasing number or Ni available sites, thus coexisting in 
both sublattices. Table IV indicates that the interaction be- 
tween Cu(Ni) and Cu(Al) atoms involves changes in energy 
of similar or smaller magnitude than the ones described in 
Section 1 . As a consequence, clustering of Cu atoms (linking 
of Cu(Al) and Cu(Ni) atoms, as in Figures 9(c), 10(c) and 
(d), and 11(d) and (e)) is observed. The NN bond [Cu(Ni) 
+ Cu(Al)]! introduces a 0.00163 eV/atom energy gain over 
the two atoms in solution, translating into Cu clustering in 
Ni-rich alloys with high Cu concentration. Finally, the small- 
est gains are realized by Ti(Al) atoms locating themselves 
at NNN distance (0.00103 eV/atom). Although small, it 
suffices to explain the observed trends in all alloys with 
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nearly the same number of Ni and A1 sites where equal 
partitioning of Ti atoms in solution in one sublattice (Al) 
and Cu atoms in the other (Ni), without interaction between 
them, is seen. 

C. Summary of Site Occupancy Behavior 

Having discussed individual and collective behavior in 
the previous sections, it is useful to summarize the main 
features by a quick review of Figure 15. While the top 
row of Figure 12 highlights the role of the coupling series 
Cu(Al) <■+ Ti(Al) ~ Ni(Al) in Ni^yjAl^-.riTi^Cu^ (* = 1, 

. . . , 5) alloys, the next lines in Figure 12 are examples of 
the same coupling feature but without the Ni antistructure 
atom, indicating that it is precisely the presence of antisite 
defects and other substitutional atoms (such as Cu or Hf [18] ) 
that ties substitutional Ti atoms together. Their absence, as 
in alloys with N A] = 33 and N M = 34, results in Ti(Al) 
atoms remaining in solution. The set with N N j = 35 marks 
the transition between the Cu(Al) +* Ti(Al) ** Ni(Al) regime, 
characteristic of Ni-rich alloys, to the migration of Cu atoms 
to the Ni sublattice that is dominant in Al-rich alloys. For 
example, Nip^Alp^TipjCufi] represents the reduced impor- 
tance of Cu-Ti coupling once Cu occupies Ni sites, while 
alloys with the same number of Ni atoms but increasing Cu 
content highlight the dual role of Cu linking Ti atoms when 
residing in Al sites and leading to Cu clustering when occu- 
pying sites in either sublattice. The NN bonds Cu(Ni)-Cu(Al) 
compete with Cu(Al)*+Ti(Al), thus explaining the small 
energy advantage of Ni [3 5 ] Al [ 33 ] Ti [ 2 ]Cu [2 ] (as shown in Figure 
12) over a similar distribution where the isolated Ti atom 
is linked to both Cu(Al) and Cu(Ni), or, similarly, the advan- 
tage of Ni [3 5 ] Al [ 32 ]Ti [ 2 ]Cu [3] over a similar configuration 
where the isolated Ti atom closes the square of NNN bonds 
between Ti(Al) and the two Cu(Al) atoms. Finally, alloys 
with N Ni = 34 and N Ni = 33 complete the transition to a 
regime where dominance of coupling between Cu atoms is 
the main feature. 

V. CONCLUSIONS 

The change in site preference in Ni-Al-Ti-Cu alloys was 
determined experimentally via ALCHEMI and theoretically 
using the BFS method for alloys. In complete agreement 
with experiment, the BFS calculations also provide an expla- 
nation for the observed site distribution of Ti and Cu atoms 
in the NiAl matrix. A simple approach, based on the defini- 
tion of local environments, provides additional information 
regarding the different features in the bonding scheme. It is 
found that subtle interactions between Ti and Cu atoms, 
together with their response to the ratio of Ni and Al atoms 
in the alloy, combine to give Ti the leading role in conserving 
its behavior when it is the only alloying addition (occupying 
sites in the Al sublattice) and give Cu a secondary role in 
individual site preference (Al sites in Ni-rich and Ni sites 
in Al-rich alloys) as a response to subtle interactions with 
Ti atoms. 


In the range of concentrations studied in the NiAlTiCu 
system, the behavior extrapolated from the ternary cases 
(NiAlTi and NiAlCu) coincides with the one observed and 
calculated for the quaternary alloys. However, additional 
features regarding the interaction between alloying additions 
can only be described by a full analysis of the quaternary 
system, which, in general, are not necessarily a direct combi- 
nation of the trends observed in the corresponding ternary 
systems. The formalism introduced in this work thus pro- 
vides an efficient tool to properly describe complex systems, 
which could be of extreme use should the interaction be- 
tween alloying additions result in features not easily predict- 
able from limited experimental evidence. 
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